Weighted Particle Variance Reduction of Direct Simulation Monte Carlo for the Bhatnagar - Gross - Krook Collision

نویسندگان

  • Donald Landon
  • David Hardt
  • Colin Donald Landon
چکیده

Direct Simulation Monte Carlo (DSMC)-the prevalent stochastic particle method for high-speed rarefied gas flows-simulates the Boltzmann equation using distributions of representative particles. Although very efficient in producing samples of the distribution function, the slow convergence associated with statistical sampling makes DSMC simulation of low-signal situations problematic. In this thesis, we present a control-variate-based approach to obtain a variancereduced DSMC method that dramatically enhances statistical convergence for lowsignal problems. Here we focus on the Bhatnagar-Gross-Krook (BGK) approximation, which as we show, exhibits special stability properties. The BGK collision operator, an approximation common in a variety of fields involving particle mediated transport, drives the system towards a local equilibrium at a prescribed relaxation rate. Variance reduction is achieved by formulating desired (non-equilibrium) simulation results in terms of the difference between a non-equilibrium and a correlated equilibrium simulation. Subtracting the two simulations results in substantial variance reduction, because the two simulations are correlated. Correlation is achieved using likelihood weights which relate the relative probability of occurrence of an equilibrium particle compared to a non-equilibrium particle. The BGK collision operator lends itself naturally to the development of unbiased, stable weight evaluation rules. Our variance-reduced solutions are compared with good agreement to simple analytical solutions, and to solutions obtained using a variance-reduced BGK based particle method that does not resemble DSMC as strongly. A number of algorithmic options are explored and our final simulation method, (VR)2-BGK-DSMC, emerges as a simple and stable version of DSMC that can efficiently resolve arbitrarily low-signal flows. Thesis Supervisor: Nicolas G. Hadjiconstantinou Title: Associate Professor

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

An improved BGK-type collision-term model for Direct Simulation Monte Carlo (DSMC) methods

The BGK and BGKC (Bhatnagar, Gross, Krook, and Cercignani) collision-term models for the Boltzmann equation are compared to the more rigorous variable-hard-sphere (VHS) collision-term model of Bird. The BGKC model is based on the well-known BGK model but uses an ellipsoidal statistical (ES) distribution instead of a Maxwellian distribution so as to incorporate both viscous effects and heat cond...

متن کامل

Nonequilibrium Bhatnagar–Gross–Krook model for nitrogen shock structure

This Brief Communication is about the study of diatomic gas shock structure based on the nonequilibrium gas-kinetic Bhatnagar–Gross–Krook (BGK) model, where an additional rotational energy relaxation equation is added. The BGK Navier–Stokes flow solver is extended to solve the nonequilibrium gas kinetic model. The numerical solutions are compared with the experiment measurements and direct simu...

متن کامل

Lattice ellipsoidal statistical BGK model for thermal non- equilibrium flows

A thermal lattice Boltzmann model is constructed on the basis of the ellipsoidal statistical Bhatnagar-Gross-Krook (ES-BGK) collision operator via the Hermite moment representation. The resulting lattice ES-BGK model uses a single distribution function and features an adjustable Prandtl number. Numerical simulations show that using a moderate discrete velocity set, this model can accurately rec...

متن کامل

Numerical Simulation of the Shock-wave Structure with Different Kinetic and Continuum Models

Gas flows in the transitional regime between continuum and free-molecular flow are widely present in nature and can be often encountered in various practical applications. Noticeable attention has been recently paid to modeling that sort of flows owing to intensive development of microand nanotechnologies and design of miniature technical devices. The Navier-Stokes equations, even with the velo...

متن کامل

Multiple-temperature kinetic model for continuum and near continuum flows

A gas-kinetic model with multiple translational temperature for the continuum and near continuum flow simulations is proposed. The main purpose for this work is to derive the generalized Navier-Stokes equations with multiple temperature. It is well recognized that for increasingly rarefied flowfields, the predictions from continuum formulation, such as the Navier-Stokes equations lose accuracy....

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2011